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SPATIAL  TRIANGULA TIO N  BY  LEAST  SQUARES  ADJUSTMENT 
OF  CONDITIONED  OBSERVATIONS 


ABSTRACT 

A  rigid  method  of  least  squares  adjustment  for  a  spatial  triangu¬ 
lation  from  angle  measurements  is  outlined*  The  method  is  based  on  an 
adjustment  of  direct  observations  which  are  conditioned.  The  results  are 
derived  for  Cartesian,  rectangular  spherical,  and  rectangular  ellipsoidal 
coordinates • 
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STATEMENT  OF  PROBLEM 


The  present  report  contains  an  analysis  of  the  problem  of  determining 
the  most  probable  position  in  space  of  a  point,  given  the  azimuth  and 
elevation  angles  to  the  point  measured  with  theodolites  at  various  stations. 


INTRODUCTION 

Because  of  errors  in  the  measured  angles  it  is  unlikely  that  any  two 
lines  of  sight  will  intersect  in  space.  Even  after  adjustment  of  the 
instruments  (determination  of  the  systematic  errors)  and  after  application 
of  corrections  to  the  instrumental  readings  according  to  those  systematic 
errors,  there  will  be  no  spatial  intersections  because  of  unavoidable  ac¬ 
cidental  errors,  which  may  be  explained  as  a  composite  of  remaining  in¬ 
strumental  and  operational  errors.  Therefore,  it  is  necessary  to  determine 
a  point  in  space  as  the  most  probable  point. 

The  criterion  for  the  most  probable  point  is  given  by  the  principle  of 
a  rigid  least  squares  adjustment $  namely,  that  the  sum  of  the  squares  of 
the  corrections  which  must  be  applied  on  the  original  measurements  becomes 
a  minimum.  This  means  that  in  the  present  case  we  must  find  that  point  in 
space  for  which  the  sum  of  the  squares  of  the  angular  (azimuth  and  eleva¬ 
tion)  deviations  between  the  observed  lines  of  sight  and  the  corresponding 
most  probable  lines  of  sight  for  all  observation  stations  becomes  a  minimum. 

Especially  if  the  measurements  are  made  by  different  types  of  measuring 
instruments,  the  observed  azimuth  and  elevation  angles  are  not  necessarily 
obtained  with  the  same  accuracy  or  in  other  words  with  the  same  weight. 
Therefore,  a  rigid  adjustment  must  provide  the  possibility  of  introducing 
weighting  factors. 


VARIOUS  METHODS  OF  LEAST  SQUARES  ADJUSTMENT 

The  purpose  of  a  least  squares  adjustment  is  the  determination  of  the 
most  probable  values  for  the  desired  quantities  by  compensating  for  the 
differences  resulting  from  unavoidable  accidental  errors.  There  are  several 
types  of  least  squares  adjustments,  depending  upon  the  measuring  method  and 
the  nature  of  the  measured  quantities.  We  may  have  cases  in  which  the  de¬ 
sired  quantity  can  be  measured  either  directly  or  indirectly,  and  where  the 
observations  are  either  independent  or  dependent.  Dependent  observations 
are  those  for  which  a  mathematical  condition  exists  between  the  observa¬ 
tions.  Accordingly,  we  have  four  different  types  of  least  squares  adjust¬ 
ment  : 

1.  Adjustment  of  independent  direct  observations. 

2.  Adjustment  of  independent  indirect  observations. 

3.  Adjustment  of  direct  observations  which  are  dependent  or  conditioned. 

4.  Adjustment  of  indirect  observations  with  condition  equations. 


PROPOSED  LEAST  SQUARES  ADJUSTMENT 


The  spatial  triangulation  of  a  point  can  be  solved  either  by  intro¬ 
ducing  the  spatial  coordinates  of  the  desired  point  (x,  y,  z)  as  unknowns 
and  by  applying  an  adjustment  according  to  the  principle  of  (2),  or  by 
considering  the  most  probable  values  of  the  measured  angles  as  unknowns 
and  adjusting  these  angles  according  to  (3),  as  direct  observations,  with 
condition  equations* 

In  any  least  squares  solution  the  main  part  of  computing  is  the 
solution  of  the  normal  equations*  To  decide  upon  the  method  means  mostly 
to  select  the  method  with  the  smallest  number  of  normal  equations*  How¬ 
ever,  further  consideration  must  be  given  to  the  coefficients  in  the  ob¬ 
servation  equations,  especially  with  regard  to  the  required  number  of 
digits*  In  the  adjustment  of  independent  indirect  observations  the  number 
of  normal  equations  equals  the  number  of  unknowns.  In  our  case,  with  the 
unknowns  (x,  y,  z),  we  would  have  three  normal  equations*  The  observation 
equations  will  not  be  linear  initially,  but  must  be  linearized  by  means 
of  the  Taylor  Series  *  This  step  calls  for  approximation  values  of  the  un¬ 
knowns*  The  required  accuracy  of  the  coefficients  can  be  expected  to  be 
relatively  high,  because  angular  deviations  cause  different  displacements 
on  the  coordinates  depending  on  the  geometrical  configuration  of  the  tri¬ 
angulation  case. 

In  the  adjustment  of  conditioned  observations  the  number  of  normal 
equations  equals  the  number  of  condition  equations.  The  number  of  con¬ 
dition  equations  equals  the  number  of  observations  which  are  in  excess . 

Because  three  angles  are  necessary  to  provide  a  rigorous  solution, 
we  have  for: 


Number  of  Stations 

2 

3 

4 

5 

n 

Number  of  Observations 

4 

6 

8 

10 

2n 

Number  of  Angles  Necessary 
for  a  Rigid  Solution 

3 

3 

3 

3 

3 

Number  of  Normal  Equations 

1 

3 

5 

7 

2n-3 

Because  of  the  fact  that  the  angular  corrections  are  small  and  of 
the  same  order  of  magnitude,  independent  of  the  absolute  value  of  the 
measured  angle,  the  coefficients  in  the  observation  equations  need  carry 
only  a  few  digits. 

The  selection  of  the  most  economical  least  squares  adjustment  may  be 
based  on  the  above-mentioned  facts* 

It  appears  that  because  of  the  simplicity  of  the  coefficients  in  the 
observation  equations,  the  adjustment  with  conditioned  observations  will 
be  advantageous  for  three,  and  possibly  four,  measuring  stations.  For 
five  or  more  stations  the  method  of  independent  indirect  observations 


* 


t 


i 


* 


l 


6 


should  be  applied.  In  missile  trajectory  measurement  the  majority  of 
missile  points  are  determined  from  two  or  three  measuring  stations.  The 
case  in  which  more  than  four  stations  are  in  such  a  favorable  geometrical 
f  pattern  that  they  will  be  used  for  a  single  intersection  is  so  infre¬ 

quent  that  it  may  be  treated  by  special  methods  such  as  grouping  of  the 
stations.  Therefore,  in  all  practical  cases  of  missile  trajectory  meas- 
\  urement  the  method  of  adjustment  of  direct  conditioned  observations  will 

be  useful. 


€ 


DERIVATION  OF  THE  SOLUTION 


1*  For  a  Cartesian  coordinate  system: 


Coordinates  of  the  camera  stations  are: 


Station 

X 

Z 

Y 

A 

Xa 

Za 

Ya 

B 

Xb 

Zb 

Yb 

C 

Xc 

Zc 

Yc 

N 

Xn 

Zn 

Yn 

The  measured  angles  at  each  station 
corrected  for  instrumental  systematic 
errors  are : 

Azimuth  :  OL  . 

Elevation:  £  . 

a 


The  elevation  angles  are  assumed  to  be  corrected  for  refraction. 


The  first-  step  is  to  establish  the  condition  equation.  For  two 
stations  with  four  observations,  we  have  4-3  =  1  condition  equation.  This 
condition  equation  can  be  obtained  if  the  elevation  Y,  of  P,  is  expressed 
first  from  a  in  terns  of  <x  ,  S  .  Y  .  <£  „  and  second,  from  b  in  terms  of 

cl  ct  ct  3. 

Sb*  ^  b* 

b  sinoc^  tan<f ,  b  sinot  tan^. 

Y  -  ha  +  ya  Sin  (c*a  +  7 a  ~  ^  +  7b  “  sin  (^a  +-  b 


h  =  3  tan  € 
a  a  a 


—  o. 


tan^., 


b  sincX^ 
sin  (aa  +«b) 


b  sincx 

a 

sin  («a  +«bj 
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and  the  rigorous  condition  equation: 


sin  o(  ^  tan  <£: a  -  sin  (X  a  tan <£  K  +  a  sin  (cXa  +  o<  K)  =  0 

Ay 


(1) 


where  a  = 


Ay  =  ya  -  yb 


and  b  =  base 


If  the  observed  angles  are  inserted  in  the  condition  equation  (1),  the 
equation  will  not  reduce  to  zero  but  leave  a  discrepancy  d. 


If  »  YoC  >  *£ 

aba 


and  V'  denote  the  corrections  which  will 


reduce  the  measured  angles  to  the  most  probable  values ,  the  condition 
equation  becomes : 

F  =  sin  (ocb  +  v^  )  tan  (<£  a  +  v£  )  -  sin  (<Xa  +  v^  )  tan  (e  b  +  v^  )  + 
b  a  a  cb 


asin(ca+v^+«b  +  v<Xb)=0 


(2) 


To  linearize  this  equation  we  have,  by  Taylor’s  Series,  neglecting  the 
terms, of  second  and  higher  power: 


d  +  (-il) 
a* 


<x. 


P 


a  f 
Sol 


Oi, 


3F 


BF 


P 


Be 


P 


Be  P 

b 


=  0 


or: 

a„ 


+  a2  v<*.  +  a3  V  +a4  V„ 
a  b  a  b 


+  d1  —  0 


(3) 


where  Y  —  180  -  (oC  y3  ) 


a^  =  -  (cos  oi  a  tan  £  b  +  a  cos  Y ) 

—  +  (cos  o(  b  tan  £  -  a  cos  Y ) 
2 

a_  =  +  sin  sec  € 

3  b  a 

2 

a,  =  -  sin  OC  sec 

4  a  b 


Li  =/0  •  d  =  /O  (sin  0<  b  tan  e  a  -  sin  c*  a  tan  b  +  a  sin  T"*) 


The  least  squares  computation  requires  now  determining  the  corrections 
v^  ,  v^  ,  v^  ,  ,  so  that  the  condition  equation  C3)  is  satisfied  and 

a  b  a  b 

that,  for  observations  with  equal  weight,  [w]  or  for  observations  of  dif¬ 
ferent  weight,  [pw]  becomes  a  minimum . 
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Introducing  a  new  parameter  A,  a  new  function  f  can  be  derived  from 
(3)  as  follows: 


f  =  P 


c* 


oc 


+  p<*  •  +  p 


oc. 


2  2 
+  P£  * 

a  eb  fcb 


+  x 


+  a2  Tc/  +  a3  +  a4  V  v  +  d>> 

3.  D  2L  D 


or,  in  another  arrangement : 

f  =  +  p «  t«2+  aix 

a  a  a 

+  P  o(  v  a  2  +  vc*  a2  * 

b  b  b 


+  p  ^  v  _  +  v.  a  A 

*a,  £a  tfa  3 

+  p  v  2+Vg-a  X 
*b  *b  %  4 


+  d»  X 

To  determine  the  five  unlaiowns  v 


,  ,  v  ,  v  ,  v  ^  and  A 
<*a  ^b 


we  have  the  five  equations: 


=  2  P<x  •  +ax* 

VCXQ  a  a 


sJr-£  ~  2  PoC  *  V<V  +  ao  ^ 
S  b  “b  2 


ST1*2  ’  vfi  +  a3  A 

£  a  a  ^ 

a 

Si=2.Pg  •  t  +a4X 


(4) 


S1  *  Vo<a  +  a2  Vo<b  +  a3  Va+&4  V£b+d=° 
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Introducing  X  =  -  2K 

a„ 


“a  V 


K 


QA  P 

-  a3 


P«. 

.  _*/ 


K 


K 


K 


C5) 


By  substitution  of  these  values  in  the  condition  equation  (3)  we  obtain 
the  normal  equation: 


[?1 


K  +  d*  =  0  or  k  =  - 


d» 

pq 

LpJ 


With  the  help  of  the  correlate  K,  the  corrections  v  can  be  calculated  from 
(5)  and  thus  the  most  probable  values  of  the  measured  angles  are  determined. 
A  partial  check  for  the  calculation  of  the  corrections  v  is  obtained  from: 


£pwj  =  -  jjd*jQ 


The  mean  error  of  an  observation  of  unit  weight  is 

f^o  —  v  num 


_  LfwI  _ 

number  of  observations  in  excess 


and  the  mean  error  of  an  observation  of  weight  p.  is  .  =  + 


A< 


/p7 


If  there  are  more  than  two  stations,  each  new  station  may  be  combined 
with  each  of  the  old  ones,  thus  giving  rise  to  additional  condition  equa¬ 
tions.  As  an  example  for  three  stations  A,  B,  C  with  six  independent  ob¬ 
servations,  we  have  6-3  =  3  condition  equations  which  are,  according  to  (3): 


(3*) 


+ 

a„ 

+  a. 

V  ^ 

+  a_ 

+ 

d« 

*a 

2 

<*b 

4 

5 

*b 

+ 

b„ 

V 

■4*  b_ 

V  , 

+  b. 

v 

+ 

d> 

°<b 

3 

OC 

c 

5 

6 

*c 

4- 

V  , 

+  c . 

v  ^ 

+  C/ 

+ 

d» 

<*a 

3 

OC 

c 

4 

6 

a 

6 

*c 

/ 


* 
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Equation  (4)  then  becomes: 

2p*_V.  +  ai  xi  +  ciA3=0 


a  a 


2ptfb  vab  +  a2  ^l  +  b2^2-0 

2V  *ck  +  b3  ^2  +  c3  ^3  =  0 
c  c 

2P£  v£  +  a4  ^  +  <>4  =  0 

a  a 

2P^  +  a_  X-.  +  =  0 

*b  b  5  1  52 

2j>€  r£  +  b6  +  c6  ^3  =  0 

c  c 

and,  for  ^  =  -  2^  /^  =  -  2K2 


3  =  -  2K3 


(4») 


a,  K3 

v: — 


a2  h  *  b2  *2 


0(,  ~  p 


'  b3  ^2  "*"  c3  ^3 

<*  =  p 

e 


a,  ^  +c,  K, 


 4 


'4  3 


*a  p* 


<5«) 


 5 


a„  Kj_  +  b„  K„ 


5  “2 


_  b6  K2  4  c6  K3 

r*c  v  v:  '  ■ 


n 


Substituting  these  values  in  (3*)  we  obtain  the  normal  equations: 


[7] 

h+[f] 

1  K2+! 

[t] 

[f] 

1  ^[t] 

K2  + 

[7] 

[7] 

1  %  +[t] 

1  M 

>] 

By  means  of  the  derived  correlates,  K^,  K^,  and  (5*)>  the  corrections 
v  may  then  be  determined 

Numerical  example  for  two  stations  A  and  B 


Measured  angles:  a.  -  38  24' 10” 

0^b  =  38  25  50 

€  =  9  06  00 
a 

=  9  43  50 
Y=cx  +a,  =  76  50  00 

a  D 


Length  of  the  base  AB  =  54  614*89 . 
Difference  in  elevation  of  the 
two  stations, 

Ay  =  +  393*80 


_Ay  _ 
~H5" 


a  =  — =  0.00721049 


Hie  elevation  angles  have  been  corrected  for  refraction: 


sin  o(  =  0*62186 
a 

sin  oi  =  0*621566 
b 

sin  Y  -  0.973712 


cos  =  0.78 
a 

cos  OC  ,  =  0.78 
b 
V- 

COS’  Y  '  =  0.23 


tan  €.  =  0.160174 

a 

tang,  =  0.171482 
b 


COS  £  =  “0.99 

a 

cos  £  ,  —  0.99 

b 


cos  £  =  0.97 

a 

cos^  £  b  =  0.97 


•—  "*  0ol3 
a^  =  0ol2 

a^  “*  +  0«6>4 
a4  “  0*64 


jaa]  =  0.85 
d!  =  +  0.000057  •  /O  =  +  11.8 


K  =  -  13.9 
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Check: 


✓ 

\ 


/ 


V 


i 


v^  =  +  1.8" 

v  =  -  8.9 

6a 

v.  =  +  8.9 

*b 


Result ; 


+  v  =  38°24>12« 

a.  O** 

cl 

c*h  +  v^  =  38  25  48 

£  +  ▼  =  9  05  51 

a  *a 

+  v  =  9  43  59 
eb 


[vvj  =  +  164.6 
-  [KdJ  =  +  164.0 


It  is  noted  that  all  computations  for  the  coefficients  a_  to  a  carry  only 

X  I4f 

two  digits.  Special  cases  of  the  solution  occur  if  not  all  angles 'have  been 
measured.  * 


If  we  have,  for  example,  three  stations,  it  may  happen  that  all  three 
azimuth  angles  are  available  but  only  one  elevation  angle.  With  these  four 
measured  quantities  we  obtain  4-3  =  1  condition  equation  and  therefore  one 
normal  equation. 

The  condition  equation  can  be  found  by  means  of  the  following  computa¬ 
tions  : 


S  sin  o L 
S  =  f1 

b  sin  (^>2  -/3) 

e  _Sbsin/3 

c  sin  -r~) 

S  sin  7* 

S  ss  c 

a  sin  or 

s±nois±n/3 sin  ir 

sin  sin  TfrW  sin  ^_1 


Introducing  od 1  =  OC+  v^ 


(6) 


=/3+v/3 
Tx  =  T  +  vr 

by  means  of  Taylor's  Series,  neglecting  the  terms  of  second  and  higher  order, 
we  obtain: 


+  cosasiiydsin  7  +  cos  ~<X)  sin  (^>2  -/3)  sin  -7“) 

+  sinc^cos^SsinT'  +  sin  (4>x  -  d)  cos  C^2  -A>  sin  (5^3  -n 

+  sin^siiy^cos  T"  +  sin  sin  cos  c«*3  -r) 

+  sinocsin^sin  T'*  -  sin  sin  (^2  -/9)  sin  c^3  -n  =  0 

ora,lV  +a,2^3  +a,3  V  +  D  =  0 

13 


C7) 


To  complete  the  analysis  the  case  should  be  considered  where  all 
three  elevation  angles  were  observed  but  only  one  azimuth  angle.  This 
case  cannot  be  expressed  in  a  simple  condition  equation.  The  practical 
computation  of  the  coordinates  of  the  missile  position  requires  the  use 
of  approximate  values  of  the  unknown  coordinates,  thus  providing  the 
means  for  an  adjustment  of  independent  indirect  observations. 


2.  Solution  for  a  spherical  rectangular  coordinate  system. 


The  measuring  stations  A  and  B  have 
the  spherical  rectangular  coordinates: 


x  ,  z  ,  y 
a  a 3  * a 


and  xb,  zb,  yb 

The  measured  angles  are: 
Station  Azimuth  Elevation 


A 

B 


oi 


The  elevation  angles  are  assumed  to 
have  been  corrected  for  refraction. 
The  length  of  the  baseline  AB  is 
denoted  by  b. 


^  \  J-  l  4-nv,  —  ,  j._  2 


r)3  (|s  +  1  tan2£:a  +  tan4*=  a^j 


>Sb[ 


Y  =  V'a2sb  tan^b+a2sb|  IT  (  °*5  +  tan2^b)+  (-1)2  tan^b  (|+tan2€b)  + 


(h + 1 tan2*  b +  tan4fi  b’l 


i,  =  (1  +  — ) 
1  N  r  J 


a2  ■  f1  + 


Ay  =  ya  -  yb 


a  - 
a3  b 


/0=  206  26498062 
alsa  tan*a  -  Vb  tan<rb  +Ay  +  alsa  M  "  Vb  Cn3  =  0 


C8) 


/ 


\ 


I 

f 
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4  • 


V 


4 


i 

1 


The  spherical  excess  of  the  triangle  A'B'P'  is  A  =  /0=  f  (cxa,cxb,  b) 


A  _  b2/0  sin*a  sinab  _  sin*a  sin0S 
“T  ^2"  *  sin  (*  a+ocb)  * 


sin  ToC+oQ 


The  values  for  can  be  tabulated  for  a  given  baseline  with  the  two 
arguments  ctf  and  Oi  ,: 


a'a=“a-T 


<*»  -o c.  - 


sin  <X! 


s.  = 


a  sxn  r'  * 


b  u 

—  -b 


b  ~T 

Yx  ~=  c<i  +a«  =oc  +a,  - 

“  '  a  h>  K 


(Set  of  Legendre) 


2A 

b  ~3~ 


sin  oC  * 


s,  = 


b  sxn  y~ 


A-b 


Substituting  s  and  s,  in  (8)  gives: 

a  D 


a1  sinO^'b  tan€a-a2  sinoc*a  tanefe+a^  sinr1  +  &1  sina'b[f}-a2  sin<X»2[lf} 

(9) 

The  most  probable  values  of  the  measured  angles  are :  (<^  )  =<X  +  v_. 

a.  a  cx 

a 

(«b}  =<xb  +  V, 


a)  =  Sa  +  T€a 

=eb  +  VS, 


(10) 


Inserting  the  measured  angles  in  the  condition  equation  (9)  we  obtain  a  dif¬ 
ference  d,  caused  by  accidental  errors  in  the  measured  angles. 

F=a1  sino£’b  (tan^a+ £l]  )  “  a2  sinc(«a  (tan<£b+[lf|  )  +  a^  sin^’  =  d  (11) 

Using  Taylor's  Series,  and  neglecting  terms  of  second  and  higher  order,  we 
obtain  from  (9)  and  (10): 


&  •&+  -~-£-  v  +  v_,  +  — v  +  —  -  v 

(  OL  a€  £a  e. 


^a  ®0(  '  b 

a  b 


=  0 


or: 


A1  v<x  '  +  A2  v«„  +  Ve  +A4v6k  +  d'  =0 

a  d  a  b 
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Hie  determination  of  ,  va  ,  v  ,  and  v^.  follows  the  procedure  out- 

a  b  6a  b 

lined  in  the  solution  for  a  Cartesian  system,  and  gives: 

d 


K  = 


M 

Check:  [w]  =  -  fKd] 


T«  =  Al‘K 
a 

v< X  ~  ^2  *  ^ 
b 

v  .  =  Aq  •  K 

^  a  3 

v  =  A,  •  K 
b  4 


A1  =  ~3oT  =  "  a2  C°Sa,a  (tan£b  +  M  +  a3  cos  r' 


-  a.  since*,  cot  Y'1 

1  D 


-jr  G-j+tan  <^a)+2  0—)^  tan€&  c|-+tan^£a)  + 


3  C'F')3  5  tan2^  a+  tan4^  a^j 

[S  1  A  S  •  8.  A  A 

~  (^+tan  £b)  +  2  -|-£tan«b  (|+ tan  <£b)  + 

2 

,  sasb  ,5.  ^  7  +  2 .  .  4  71 

3  “T“  5  tan  £  b+  tan  €  b^ 


A2  =  -fi  =+al  cos0<r,b  (tanea  +  [i]  )  +  a  cos^* 


+  ab  sinoc'b  cosec  r»  (^+tan2£a)+2  -|Jitan£a  (|+tan2€a)  + 

2 

s  *  s 


s»s, 
a  b 


A  -  3F  -+ 


,  a  b  ,5  ^  7  ..  2^  ,  .  4  ,"I 

3  -j-  (jj+  E  tan  <s a  +  tanV^J 

[s  s 

(^+ tan2<?b)  +  2  (^)2  tan«b  (^+tan2^b)  + 

3  (-F)3  tan2€  b+tan4eb^ 

n  sino(*  T"  s  s  0 

i-,— S  |l+  (-|)  2  tan*a+  (^)2  (|+3  tan2£a)  + 

a  ^ 

(-jr)3  (|  tane&+  4  tan3^a)^J 


4 

f 
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.  _  3F  _  a2  sinot'a 

9*---zzrr 


”  S  S 

1+  2  tan€b+  (^)2  (|+  3  tan2£b)  + 


(~)^  tan£‘b+  4  tan^S  b)J 


Because  of  the  fact  that  the  coefficients  A.  —  A.  must  be  computed  with  a 

-L  4 

few  digits  only,  considerable  simplifications  are  possible.  The  follow¬ 
ing  formulas  may  be  considered  adequate  for  ranges  up  to  150  miles  and 
elevation  angles  up  to  80°,  even  for  high  precision  instrumentation.  In 
the  follovdng  formulas  the  number  of  terms  to  be  carried  depends . on  the 
geometry  involved  in  the  triangulation  case  and  on  the  accuracy  of  the 
instruments  used. 

The  difference  d  is  calculated  from  (11)  as  follows: 


d  =  a^  sinoe»b 

(tantfa  +  [i]  )  - 

a,,  since1 

2  a 

(tan  £  b  +  [ 

„  _  ,  .  7a 

an  — ,  1  Hr  — 

yb 

a_  =  1  +  — - 

a  =  Ay 

II 

<<! 

P 

1 

O' 

1  r 

2  r 

3  b 

to 

<x*  =  oe  -  4_ 
a  a  3 

oe»  =a,  -  -4- 
b  b  3 

r>  =  oe  1 
a 

+  (Y  t 

b 

b  sintf’-j 

s  = - 1 — — p 

a  sin  jr « 


Sb  = 


b  sinot 
sin  r'1 


A  since  since, 

A  _  K  a _ b 

*3~  sin  (ex  +oe,  ) 

a  b 

h2 

K  =  £-£L;/0=-206  265 1 
6r 


can  be  tabulated 


[l]=T-(|+tan2£a)+(^.)2tansa  (§+ tan2£a)+ (!i)3  (jj+  5  tan2£a+  tan4£a) 
L11]  =-r  (j+ tan2eh)+  (-i),2  taneh  (jj+  tan2«= J  +  (-1)3  (|j+  5  tan2ek+  tan4fij 


A-L  =  -  a2  cosot^  (tan^b  +  jXlj  )  +  a^  cost'* 
A2  =  +  a^  cosce'b  (tan<£a  +  QlQ  )  +  a^  cos  T-*1 

S 


a..  sin  ex’,  s 

A0  =  +  - -  (1  +  2  tan<^a) 


- =ZT 

cos  <£ 


a^  cos  ?-*  is  only  sig¬ 
nificant  if  the  differ¬ 
ence  in  elevation  be¬ 
tween  the  measuring 
stations  is  considerable 


a  sinoe*  s. 

A  ,  =  £ -  (1  +  2  —  tan^O 

cos  ^ 
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In  case  all  the  elevation  angles  are  not  available  the  condition 
equation  (7)  may  be  used.  The  angles  must  be  changed  first  to  plane 
angles  by  subtracting  one  third  of  the  corresponding  spherical  excess . 

3®  Solution  for  an  ellipsoidal  rectangular  coordinate  system 

The  radius  of  curvature  for  a  point  on  the  ellipsoidal  earth  changes 
with  the  azimuth.  If  M  is  the  radius  in  the  meridian  and  N  the  radius 

perpendicular  to  the  meridian,  there  is  the  relation  ^  =  Y  .  The  more 

nearly  this  ratio  is  equal  to  one,  the  more  nearly  may  the  earth  be  con¬ 
sidered  as  a  sphere  in  the  vicinity  of  the  point . 


2 

The  values  V  are  rather  different  from  1 

o  2 

and  only  on  the  poles  &  =  90  )  does  V 
become  equal  to  1.  For  a  latitude  of  30° 

7/e  have  ,  a  deviation  of  about  l/2$. 

Because  Tire  may  neglect  for  our -purpose  the  change  of  curvature  with' 
the  change  of  the  mean  latitude  of  the  triangulation  side  involved,  it 
is  possible  to  substitute  for  the  ellipsoid  at  a  measuring  station  a 
number  of  spheres  whose  radii  change  with  the  azimuth  of  the  considered 
sides.  The  radii  can  be  tabulated  for  a  station  at  suitable  intervals. 
The  ellipsoidal  calculation  may  now  use  the  same  formulas  which  were 
derived  in  the  chapter  on  a  spherical  coordinate  system.  The  only  dif¬ 
ference  is  that  instead  of  a  single  radius  for  all  measuring  stations 
a  specific  radius  for  each  line  of  sight  must  be  introduced. 


<P 

Latitude 


S-v2 

i 


30; 

45 

60c 

90c 


1,007 

1,005 

1,003 

1,002 

1,000 


NUMERICAL  EXAMPLE  FOR  THREE  STATIONS  IN  AN  ELLIPSOIDAL  RECTANGULAR 
COORDINATE  SYSTEM  'WITH  DIFFERENT  TYPES  OF  INSTRUMENTS  AND  INCOMPLETE 

ANGLE  MEASUREMENTS 


The 

ellipsoidal  rectangular  coordinates  of  the 

measuring  stations. 

expressed  in  meters. 

are: 

% 

Station 

x  . "  z 

7 

N 

+  37010.871  -  31383.845 

1524.003 

P 

-  11088.230  0 

1219.202 

0 

+  92519.561  +  46815.110 

2133.604 

The 

length  of  the  baseline  NP  is  b  =  57432.138 

Ay  =  yn  -  yp  =  +  304.801 

a3  =  TT  =  +»°°5307 
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\ 


i 
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After  the  measured  angles  are  corrected  for  systematic  instrumental 
errors  according  to  adjustment  charts  and  the  refraction  correction  is 
applied  to  the  elevation  angles,  we  have  the  '’measured”  angles: 


Elevation 

Station  Horizontal  Angle  Angle 


Instrumental 
Type  of  Accuracy 

Instrument  .  (Mean  Error) 


N 

P 

0 


112  59' 45"  from  MP 
31°00’00”  from  PN 
20  05 '  —  from  OP 


80° 00 ’ 10”  Cine-t  heodolite 

72°29'.54.0«  Ballistic  Camera 
— ..  .  Tracking  Telescope 


+  20” 

+  3" 

+  200” 


If  we  assign  unit  weight  to  a  Cine-theodolite  measurement,  the  following 
weights  must  be  introduced: 


Station 

N 

P 

0 


Weight 


Because  of  the  fact  that  we  have  no  elevation  information  at  station  0, 
we  have  five  measured  quantities  or  5-3  =  2  condition  equations* 


Erom  the  coordinates  of  the 
stations  the  angles  ^ 

are  known*  They  have  already  been 
corrected  for  spherical  excess. 


=  92014'39*1" 
^>2  =  30°18’56*7” 
=  57°26’24.2” 

£=180o00’00.0” 


The  condition  equations  are,  from  (11)  and  (7): 

/°[ai  sintf'p  (tan^n+ [ij  )-ag  sino('n  (tang  )  +  a^  sin 7^’  =  d» 

/ofsin  {<£>  -&.'  )  sin  +OL'  )  sinoc'  +sin(X’  sinoi'  sin  (<^>  )  = 

'  L  W1  n  r  2  o  p  n  o  3  p 


al=1+r 


n  _ 


n 


1.00  240  r  =  6363154.8  (for  azimuth  (NP)  -0 (  ) 
n  n 


D» 
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y„ 

a2=l+-ii=  1.000  192 

rp 


r  =  6353492.0  (for  azimuth  (PN)  +  a  ) 

p  p 


a3 

_ 

“  IT 

—  +0© 

005  307  The 

values 

n 

Si 

H 

A  _ 

“T  " 

112°59,42.7“ 

A 

“r 

p 

II 

R 

A  _ 

3 

30°59!57.7” 

■V 

r* 

=  00 

n 

+a*  = 

p 

:  143°59!40.5" 

00 

n 

=  20°05 

1 

The  values  for  r  can  be  taken  from  a  table. 


=  2.3"  can  be  taken  from  a  table 

■with  the  arguments  Oi  and  oc  . 

n  P 


The  spherical  excess  for  oC 

can  be  neglected  due  to  the 
mean  error  of  +  200". 


00 


n 


sin  +  0.920  5375 

COS  -  0.391 


p 

+  0.515  0287 
+  0.857 


+  0.587  8618 
-  0.809 


„  sm 
3  *  cos 

+  0.003120 

-  0.004 


r« 


b  sin  OC  * 


sn  = 


s  = 
P 


sin 

b  sin  04* 

n 

sin  t"1 . 


P  


50  316.587 


=  89  933.445 


(Set  of  Legendre) 


(-»)  = 
vr  y 

n 

0.007  907 

tan  € 

S  o 

C— )  = 

Vr  J 

n 

0.000  063 

tan2£ 

11 

„nUH 

0.000  0005 

tan*£ 

H  - 

0.270  667 

tan  £ 

= 

0.014  155 

tan 

P 


=  5-672890 

=  32.181682 

=  1035 .6607 

+  [l]=  5.943  558 
=  3.171  273 


i  +  tan2«=  =  32.681  682 
<.  n 


f.  +  tan2e  =  33.015  015 
o  n 


+  1  tan2€  n  =  37.753  627 


i  +  tan2  =  10.556  973 

£  P 


V 


i 

( 
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»_  5 

(-£)  =  0.000  200 
rp 


tan2^  =  10.056  973  |  +  tan2<£  =  10.890  306 

p  bn 


(-H)-3  =  0.000  003 


tan4£  =  101.1427 
P 


^  +  |  tan2*  p  =  ^*941  469 


|jl]  =  0.156  674 


tan  €  +  |II|  =3.327  948 

P 


cose  =  0.174 
n 

COS£.  -  0.301 
P 


COS  =  .  0.030 

2 

cos  e  =  0.090 

P 


d  =  0.000  867 

d»  =  +  178.83 


A,  =  +1.30  A0  =  +  5.09  A  =  +  18.63  A  =  -  11.09 


2 

>3 


Sin 

Cos 

<Xt 

n 

=  - 

20°45,03.6“ 

-  0.354 

308 

0.935 

00 

0 

=  + 

50o23! 56.7ff 

0.770 

503 

0.537 

<X* 

p 

=  + 

26°26  •26.5s 

0.445 

271 

0.895 

00 

n 

=  112°59,42.7w 

0.920 

538 

-  0.391 

P 

= 

30°59’57.7« 

0.515 

029 

0.857 

OO 

0 

=  - 

20°  Q5 '  — 

-  0.343 

387 

0.939 

D  =  + 

.0001500  D«  =  + 

30.94 

B1 

=  +  0.43 

b2  =  +  0.27 

=  -  0, 

and  the  condition  equations  are: 


+1.3°  +  5.09  +  0.  va  +  18.63  vg  -  H.09  *  +  178.83  =  0 

p  o  n  ^ 


n 


+0..43  -  0.52  +  0.27  +0  +  0  v  ^  +  30.94  =  0 

n  p  o  n  p 


[ir]  =  351.72 


M- 


=  +  0.51 


[?]- 


+  7.48 


and  the  normal  equations: 
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351.72 


+  0.51 
+  7.48 

v  =-2^" 

n 

V  =  +  0.0 

P 

va  =  -  111. S'7 

o 

y2l_&  2 

^ =  +  10.4”, 

and  the  mean  errors  of  the  observation  are: 

for  station  N:  +  10.4" 

0:  +  104.0" 

P:  +  1.5" 

If  we  apply  the  corrections  to  the  measured  angles  we  obtain  the  adjusted 
angles: 


Kg  +  178.83  =  0  Kx  =  -  0.50 
Kg  +  30.94  =  0  Kg  =  -  4.14 


ve  =-9.3” 
n 


v^  =  +  0.1 

S 


U 


Check:  Jjpvv]  =  217.7 
f  kdj  =217.5 


an  =  112°59,42.6"  e  =  80°  00’ 00.7"  * 

<X  =  20O03^08.2» 

0 

«p  =  31°od«00.0"  =  72°29»54.1“ 

These  adjusted  angles  represent  now  a  rigid'  solution.  An  independent  check 
for  the  adjusted  angles  may  be  obtained  by  substituting  the  adjusted  values 
in  the  condition  equations  (6)  and  (9)>  which  must  be  satisfied. 


CONCLUSIONS 

The  adjustment  with  conditioned  observations  is  useful  for  a  spatial 
triangulation  from  two  to  four  measuring  stations. 

The  method  enables  us  to  use  incomplete  station  measurements  and  to 
combine  different  instrument  types  by  introducing  weighting  factors. 

The  accuracy  with  which  the  result  is  obtained  depends  on  the  number 
of  terms  carried  in  the  formulas.  The  formulas  are  expressed  in  series 
form,  thus  providing  a  means  for  breaking  off  the  numerical  work  at  the 
most  economical  point. 

The  method  is  a  rigorous  adjustment  •  The  calculated  corrections  for  ■ 
the  measured  angles  are  obtained  without  recourse  to  approximation  values 
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of  the  coordinates.  Information  about  the  consistency  of  each  of  the 
measuring  instruments  is  directly  available  by  comparison  of  the  most  , 
probable  angular  corrections.  A  calculation  check  is  provided  by 


After  applying  the  corrections  to  the  measured  angles, 


the  triangulation  presents  a  rigid  solution  and  the  coordinates  of  a 
point  may  be  calculated  from  any  two  stations.  In  this  connection  a 
proposal  for  an  ellipsoidal  rectangular  coordinate  system-  is  outlined  in 
another  report^-. 
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